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The QCD (tri)critical point, a genuine second- order phase transition, implies existence of a massless 
scalar mode, and singular behavior near it. In this work we however focus on the finite region around 
it, defined by a condition m a — 2m w . We point out that in this region the inter-hadron interaction 
should be dramatically changed. Light sigma should increase attractive mean field potentials for 
baryons and non-Goldstone mesons. The same effect can be observed in additional downward shift 
of the mass of vector mesons p, u>, <f>, accessible via dilepton experiments. For pions we predict that 
the mean change due to light sigma is in fact a repulsive mean potential. The implications of these 
effects for collective pion and nucleon flows, radial directed and elliptic, are estimated. Finally, we 
speculate tha unusual behavior of flows observed by NA49 at 40 GeV PbPb collisions at SPS may 
be explained by location of the critical point region nearby. 



I. INTRODUCTION 

A. The QCD (tri)critical point 

Let me start by briefly reminding well known facts. Ig- 
noring strangeness and starting with massless u,d quark 
theory, one arrives at the situation in which the phase 
transition at small baryon chemical potential /i is a sec- 
ond order line, changing to a first order line at the so 
called tricritical point. In this theory it is convenient to 
look at pions and sigma as one 4-component field <pi and 
write Landau-Ginzburg potential in terms of its square 

^ = \{Mi) + \{4>i<t>i? + l^i4>if (i) 

The coefficients a(T,mu),b(T,mu) 7 c(T,mu) are some 
functions. The second order line is defined by zero mass 
condition a(T, mu) = 0: at it all 4 fields tt, a are massless, 
so the transition belongs to the universality class of the 
0(4) spin model. At the tricritical point additionally the 
second equation holds b(T, mu) = which together with 
the first one select one critical point. Here the indices 
are as in the mean field theory. 

In the real world the nonzero quark mass adds term ma 
(and other odd terms) to this nice symmetric Lagrangian 
and makes it a bit more compacted. Still straightforward 
manipulations with this function tell us how all quantities 
change along the critical line, see the second paper [1]. It 
makes all fields massive in the vacuum, and transform the 
second order line to just a crossover (see Fig.l ). ri The 
tricritical point is changed to just a critical point with 
Ising exponents and just one massless mode, m a = 0. 
The pions are massive there, and m\ ~ m 4 / 5 which is 
only slightly differs from m 1 in G-Oaks-Renner relation 
in the vacuum. Below we will ignore small correction 
to it, ~ (Aqcd/itI') 1 ^ 10 , and assume the pion mass at 
the critical point is the same as in vacuum. Obviously, 
there should be a finite region (the shade area in Fig.l to 
be refered to as "oval" below) in which m a < 2m^ . One 
consequence of this condition is that a — > 2tt decay is im- 
possible: but the boundary of the oval will be important 
for other reasons as well, see section IV. 



Stephanov, Rajagopal and myself [1] were the first to 
propose experimental search for the QCD (t)critical point 
in heavy ion collisions, by varying the energy of and look- 
ing for "non-monotonic" signals. The first specific signal 
proposed in that paper was an increased event-by- event 
fluctuations near it, reminiscent of critical opalescence 
well known near other 2-nd order phase transitions. 

Although well motivated, the fluctuations is a rather 
subtle signal, because (as in fact emphasized in original 
papers on event-by-event fluctuations [2,3]) the experi- 
mentally observed fluctuations are most likely to be just 
equilibrium thermodynamical fluctuations at the ther- 
mal freeze-out. Therefore, critical fluctuations can only 
be seen due to some non-equilibrium "memory effects". 
Asakawa, Heinz and Muller [4] and Jeon and Koch [5] 
suggested that large-scale fluctuations of conserved cur- 
rents (such as electric and baryon charges) are more slow 
to dissipate, and argued that those may show the "pri- 
mordial QGP" values, rather than those for hadronic 
gas at freezeout. Unfortunately further studies (e.g. by 
Stephanov and myself [6]) have shown that the diffusion 
from the boundary of experimentally covered begins to be 
large enough to wipe it out. Experimentally, at any col- 
lision energy studied so far, the charge fluctuations seen 
are close to resonance gas values, while small expected 
"QGP values" unfortunately are not observed. 

The second proposal of [1] was the so called "focus- 
ing effects", a particular deformation of adiabatic cool- 
ing lines in the vicinity of the critical point. Nonaka and 
Asakawa [7] have provided a description of this effect in 
significant detail. 

In this paper I would like to reconsider the issue and 
suggest (hopefully) more robust signals of the critical 
point, which would not rely on subtle memory effects. 
One of those is collective flows, which show the effect ac- 
cumulated during the whole evolution along the cooling 
path, rather than just its end the freezeout. The second, 
even more direct but more challenging observable, is the 
invariant mass spectrum of dileptons, some of which pro- 
duced right when the system's cooling path passes the 
oval, perhaps even near the critical point. 
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FIG. 1. Schematic view of the critical point, for 2 fla- 
vor QCD with nonzero quark masses. The black dot mark 
the position of the critical point, separating a crossover (the 
dashed line) from the first order transition (solid line). The 
shaded oval shows the region where m CT < 2m w and its decay 
into pions is forbidden, while the wider ellipse corresponds 
to the region in which it is less than twice the average pion 
energy,m CT < 2 < E V (T) . 

The reduction in sigma mass near the whole critical 
line, and especially at the beginnig of the crossover re- 
gion \i — 0, T = T c were discussed in numerous models 
since 1960's, see e.g. [8] and references to earlier works 
therein. But susceptibilities obtained in lattice simula- 
tions before, at finite T and zero ix, do not show any 
large peaks corresponding to light sigma* . At the critical 
point, however, the sigma mass must be zero by defini- 
tion, regardless of the quark masses^. Furthermore, this 
work was triggered by new lattice data (see next sub- 
section) have now indicated a quite different behavior of 
these susceptibilities, consistent with quite small sigma 
masses and consequently large effects at /i w T c . 

The discussion of critical properties themselves reduces 
to the statement that it belongs to the simplest Ising 
universality class, with standard critical exponents. The 
issue of dynamical exponents and universality class was 
discussed by Son and Stephanov [9]: the only relevant 
conclusion the reader should be reminded now is that 
the dynamical index z, relating the available time the 
system spends near the critical point 
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with the maximal correlation length £, is z « 3. There- 
fore, in a realistic central collisions of the heaviest nuclei 
available, the time limits the correlation length and the 
sigma mass by about 



Still, a significant reduction of the sigma mass from its 
vacuum value, by about an order of magnitude, may be 
possible. As we argue in this work, even a decrease by 
about factor 2, bringing us close to "oval" in Fig.l, should 
result in quite dramatic changes in inter-particle interac- 
tion. 



B. Recent lattice results 

Discussions of the possible search for the critical point 
in experiment has resulted in a search for it on the lattice 
as well. Of course this is a very difficult task, as direct 
simulations at finite /i are impossible due to the notorious 
sign problem. Two methods pursued are re-weighting 
along the critical line [10] and Taylor expansion in powers 
fx/T. (We remind the reader that \x in this work is defined 
per quark, so fx/T = 1 means the usual baryon number 
chemical potential to be about 500-600 MeV.) 

Let me just focus on one most resent study of the lat- 
ter type, by Bielefeld- Swansee group [11] which followed 
such an expansion for the 2-flavor QCD up to terms 
0((fi/T) 6 ). The authors conclude that they don't ac- 
tually see any direct signatures of the critical point yet, 
nor in fact any other indications for the presence of the 
first order transition at larger [x. In the covered domain, 
namely lx/T < 1, the Taylor series seem to converge rea- 
sonably well. This is not surprising, since their quark 
mass is still quite large compared to physical one, which 
is suppose to move the critical point toward larger \x. 
And yet, a closer look at their results show that their 
numerical data are quite remarkable, they show surpris- 
ingly abrupt changes in the system, and some indications 
are perhaps relevant for the issues discussed in this work. 



*For the quark mass value used in the calculations, of course. 

^Note that although we know from experiment how experi- 
mental sigma resonance is seen in pion and kaon scattering, 
we do not know that for the critical mode. Thus its interac- 
tion with strange quarks may be different from that with light 
ones: the corresponding mixing angle between SU(3) singlet 
and octet is unknown parameter. 
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FIG. 2. The quark number susceptibility Xq/T 2 (left) and 
isovector susceptibility \i /T 2 (right) as functions of T/To for 
various (X q /T ranging from fi q /T = (lowest curve) rising in 
steps of 0.2 to lAq/T — 1, calculated from a Taylor series in 
6 th order. Also shown as dashed lines are results from a 4 th 
order expansion in fJ. q /T. 

In Fig. 2 borrowed from this work we show the nor- 
malized quark number susceptibility Xq/T 2 , for the 
usual (baryon number) chemical potential, as well as for 
isospin-related one. One can see that, as fi/T grows, the 
former one develops a sharp peak below T c , while the 
isovector susceptibility does not. 

Is this peak due to 4-th order 2-loop diagram with a 
sigma meson exchange, shown in Fig. 3 (a)? It would show 
a peak in /i 4 coefficient (as indeed observed on the lattice) 
and no peak in the isovector case (again, as observed)? 

Redlich and Karsch [12] argued that (at least the 
left-hand side of the peak T < T pea k) is not due to 



this diagram, which is suppressed by a need to excite 
a baryon in each loop and is thus suppressed by extra 
~ exp(-M N /T) « 1/300. They find it hard to think 
that even enhancement due to small sigma mass ~ 1/ m 2 
of the diagram (a) can beat that. They argue instead 
that the l.h.s. of the peaks is quite well described by a 
baryonic resonance gas^ and the usual one-loop diagram 
(b). 

If so, why does the signal drop so rapidly at the right- 
hand side of the peak? One possible answer is rapid melt- 
ing of many baryonic states above the critical line. An- 
other possibility is that the susceptibility is decreasing 
because (in spite of deconfinement) the baryon mass is 
larger than in vacuum. In fact, the effective quark mass 
(defined as half potential between qq at large separation) 
is indeed so large above the critical line, that quark con- 
tribution ~ exp(—Mqff(T)/T) can get even more sup- 
pressed than that of the nucleoli. 
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FIG. 3. Two diagrams responsible for 4-th order part of 
the density-density susceptibility. The black dots stand for 
the original vector density operators, the open circles are for 
insertions of the mean vector density of the matter. The solid 
lines are states with the baryon number, the baryons below 
the critical T and quarks/diquarks above it. The dashed line 
represents the a meson. 

Now let us look at the behavior of the chiral condensate 
and chiral susceptibility snown i n Fig-4 from the 

same lattice study. The condensate shows little change, 
just slight shift downward in the critical T as fi grows, 
as expected. The susceptibility however seem to have a 
larger peak: in contrast to plots shown above, now there 
is no need for extra vector current insertions, so it prob- 
ably is the (unsuppressed) sigma contribution. Thus, in 
spite of larger uncertainties and similarity to the plot 
above, I still think this peak does indeed contain a sig- 
nificant part of a sigma exchange, and therefore the peak 
is the first manifestation of reduced sigma mass. If the 
absolute hight of the peak is a measure of l/m^(/x), one 
may thing the reduction of the mass itself is by a factor 
2 or so. Thus, optimistically, one may think that at least 
the edge of the "oval" of Fig.l is more or less reached 
by those simulations. Needless to say, more work is need 
to see if this is correct: the simplest of those is to see a 
correlator rather than integrated susceptibility. 



*With masses of mesons and baryons correctly adjusted to 
quark masses used in the lattice simulations. 
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FIG. 4. The chiral condensate (left) and chiral sus- 

ceptibility Xi>^ (right) as a function of T /To for \i q jT = 0, 0.4 
and 0.8. The chiral condensate drops with increasing n q /T 
and the peak in becomes more pronounced. 



Summarizing this subsection, the results of numer- 
ical simulations presented in [11] does not found the 
critical point. Instead, they show a dramatic change 
of the baryon contributions to thermodynamics, from a 
resonance gas with unmodified masses, to a completely 
different regime in QGP. Chiral accessibility possibly 
show reduction in the sigma mass by a factor of 2 for 
fi w .8T C ~ 130 MeV, which however remains rather un- 
certain. 



in textbooks in detail. We will not discuss spin-spin or 
spin-orbit parts of it, but focus rather on one model, as 
simple as possible, for the central forces. A well known 
Walccka model [13] is sufficient to demonstrate our main 
points. 

Those forces can be viewed as a combination of the 
sigma exchange^ and the omega exchange terms 
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In this simple model g s = M n j w 10 and g\ = 190 
which fits the data well enough. 

The most important lesson about nuclear forces we 
would like to remind the reader is that the nuclear po- 
tential is in fact a highly tuned small difference of two 
large terms. Although both terms are comparable to 
the nucleon mass, the resulting potential has a min- 
imum of only about a percent of the available mass, 
V(r min )/2M N ~ O(10~ 2 ). ** The so called "relativistic 
systematics" , the dependence of nuclear forces on relative 
motion, reveals this fact well, as the scalar and vector 
parts change differently under Lorentz transformation. 
It makes the forces repulsive at semi-relativistic energies, 
as multiple NN and heavy ion low energy collisions have 
well documented. 

Because of this fine tuning, the nuclear forces are quite 
sensitive to the 4 parameters of the model, two couplings 
and two masses''"!'. 

If one simply takes m a — > 0, keeping the other param- 
eters the same, one gets huge attraction V ~ —lGeV. If 
on the other hand the omega mass is also goes to zero, 
the result is huge repulsion^ V ~ +lGeV . More mod- 
erate version of the same exercise (with parameters to be 
explained later) is demonstrated in Fig. 5 . 



§ Strong variation of the phase shift at 400-600 MeV is known 
for a very long time in tttv scattering, and was also identified 
in the attractive part of the nuclear forces, see e.g. a review 
on applications of Walecka model [13]. Whether one would 
like to call it a resonance or not, the fact remains that it 
dominates the attractive part of NN interaction, and is thus 
responsible for nuclear binding and for our very existence. 

**Of course, in nuclei there is another cancellation between 
kinetic and potential energy, leading to a typical nuclear 
physics scale to 1 MeV ~ 10~ 3 Mjv. 

^Strong sensitivity of the deuteron binding to and even- 
tually to the ratio of quark masses to Aqcd , combined with 
the Big Bang Nucleosynthesis data, makes the best limits on 
the cosmological variation of these quantities [14]. 

^Because the omega coupling is larger than the sigma one. 



II. NUCLEAR FORCES 



Studies of the NN forces, are at the very basis of nu- 
clear physics, go back by many decades and are discussed 
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FIG. 5. The NN potentials in the Walecka model, for 
the following values of the masses: m CT ,m w — 600,770 MeV 
(black), 280,770 (red), 280,500 (blue) 

A nucleon in a homogeneous matter is in a mean field 
potential which can be easily calculated in this model. 
The modification of it due to modification of the a- 
exchange diagram^ 



q* 2 Q 2 
AV = -n s {^-^) 
m* m% 



(5) 



where n s is the total scalar density of matter. Unlike 
vector density, it includes all quarks and antiquark with 
the positive sign, so 



n s = n B + n s + -n n Gr, 



(6) 



where nGm stands for all non-Goldstones mesons. The 
factor 2/3 corresponds to a simple additive approxima- 
tion, assuming that the scalar coupling gl is additive for 
valence quarks. 



III. VECTOR MESONS 

The mass shifts for p, u due to sigma exchanges are of 
the same nature as for baryons, so in the same spirit of 
additive coupling one expects extra*** mass shift to be 
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§§ What exactly is a sigma field and how it is interaction with 
pions and nucleons is described e.g. in [15] 

***Which obviously goes on top of other effects, such as res- 
onances and the sigma exchanges with unmodified mass. 



These shifts were experimentally seen in heavy ion col- 
lisions e.g. in peripheral AuAu collisions at RHIC, by 
STAR collaboration [23] for few resonances, including 
vectors p and K* . In the former case the mass shifts 
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-70 MeV, with about half of that esti- 



mated to come from the sigma exchange [22]. If so, one 
may thing that at the edge of the "oval" this term would 
be 4 times larger, or about — 120 MeV. Very close to 
the critical point the effect is still finite, in spite of zero 
sigma mass, due to vanishing coupling: this however only 
happens when sigma size gets really large. 

The resonances like p seen via irir channel arc of course 
observed near the kinetic freezeout conditions, that is at 
rather dilute matter. One may wander if those are not 
too far from the critical line and critical point, to show 
any observable effect. 

Another well known possibility is observation of vector 
states via dileptons: in this case they are collected not 
from freezeout but from all the 4-volumc of expanding 
fireball, with the inside of the "oval" included. Since we 
expect additional shift of vector mass, the enhancement 
is supposed to increase roughly as exp(AM(T)/T) <~ 3 
due to twice lighter sigma. 

We do indeed know that NA45 (CERES) collabora- 
tion have seen rather large enhancement of small mass 
dileptons at SPS, especially at 40 GeV. Recent run of 
NA60 experiment with dimuons is expected to provide 
additional clarification on this issue. A non-monotonous 
dependence of the light dilepton continuum on collision 
energy would be possible indication for the critical point. 



IV. THE INTERACTION OF THE PIONS 

Unlike baryons and non- Goldstone modes of the the- 
ory, the pions cannot simply linearly interact with the 
sigma field, because any nonzero scalar density would 
made them massive even in the chiral limit and thus vio- 
late the Goldstone theorem. In the original sigma model 
formulation that is resolved by some cancellation between 
diagrams. Better solution is to change the scalar field to 
its "radial" version on the chiral circle (see [15]), after 
which sigma interacts only with a derivatives of the pion 
field and preserves the theorem in each vertex. 

As in the chiral limit the interaction of low energy pions 
is described by Weinberg Lagrangian, 
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the corrections due to nonzero scalar density of the 
medium can simply be absorbed into a modified coupling 
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to be substituted into known results about pion gas ther- 
modynamics [17,18] and kinetics [19]. For low-T pion gas 
the same is true for realistic QCD with massive quarks 
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and pions, since 0(m) corrections to scattering (such as 
Weinberg scattering lengths) are very small. 

A real hadronic matter produced in heavy ion collisions 
have sufficiently high temperature T > Tf <~ 100 MeV to 
excite resonances, which truly dominate the interaction. 
For pions those are a and p. those provide scattering 
rate and ImV. The real part of the mean field potential 
for pions is very small and usually neglected. 

The (momentum-dependent) potentials are usually de- 
fined as 

V(p) = (m 2 + p 2 + £(;p)) 1/2 - (m 2 +p 2 ) 1 / 2 (10) 
where 

^£/|^^> (") 

i 

is the mass operator related to the forward scattering 
amplitude M. The sum stands for all particle types in 
matter. 

The amplitude can in turn be written as a sum over 
resonances in each channel, and we are going to focus on 
a in 7T7T. For this case the forward scattering amplitude 
is of standard Brcit-Wigner type 
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where E = ^fs and q are the CM energy and momentum 
of two colliding pions with momenta p, k. We used the 
total width in the form T tot = (AGeV)(q/E) where q, E 
are the CM. momentum and energy per pion: not that 
it vanishes at q=0, as is needed by the phase space. The 
Tin = rt t/2. 

As it is well known, the real part of the amplitude 
changes sign at the resonance, and the sign means that 
effective potential obtains positive ( repulsive) contribu- 
tion from 7T7T states above the resonance, E > m CT , and 
negative (attractive) from the states below it. These two 
parts of the integral over k tend to cancel each other. As 
a result, the interaction between pions is weakly attrac- 
tive in a pion gas, with a potential strongly dependent 
on the pion momentum p. 
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FIG. 6. Effective potential for a pion at 
rest Re[V e ff(p = 0)] [GeV] induced by sigma resonance, as 
a function of the sigma mass M a [GeV] . 

Another (somewhat more intuitive)) picture of pion 
in a pion gas is that since they spend part of the time 
rotating around each other in form of resonances, the 
average pion's velocity is slightly reduced in the pion gas 
relative to that in the vacuum [20]. 

However this is going to change near the QCD criti- 
cal point. Near the dashed line in Fig. 1 the effective 
pion potential rapidly changes and is getting repulsive 
inside the solid shaded oval, as seen from Fig. 6. For 
light enough sigma, all colliding pions are always above 
the sigma resonance, and thus the attractive contribu- 
tion is absent, with only repulsive contribution of sigma 
resonance remaining. As seen from the figure, it leads 
to a rapid change of the potential, which is not small. 
It happens exactly at the boundary of the oval, where 
m a = Im-x and the scattering length of two pions at rest 
goes to infinity since the amplitude is ~ 1/q and the rela- 
tive momentum q — > here. The corresponding tttt cross 
section at this point is reaching its s-wave unitary limit 
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This behavior is identical to the so called "Feshbach res- 
onance" at zero energy, which is used in atomic physics 
for cold trapped atoms. The result is known to be a new 
"strong coupling" regime of matter in which it shows 
liquid-like behavior [21]. 

In summary of this subsection: we predict that close 
to QCD critical point the pion-pion interactions rapidly 
change their sign, leading to repulsive mean potential. 
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FIG. 7. Elliptic and directed flows of pions and protons 
versus rapidity at 158 A-GeV Pb+Pb collisions [16] measured 
for three centrality bins: central (dots), mid-central (squares) 
and peripheral (triangles). The solid lines are polynomial fits 
to the data [16]. 



V. DIRECTED FLOWS 

One general difference in the vicinity of the critical 
point is that near it the hydrodynamical description is 
to be appended by simultaneous solution for the critical 
slow mode, the a field, coupled to matter evolution***. 
The first pioneering attempt to do so has been done re- 
cently by Dumitru, Paech and Stocker [26]. 

We have argued above that near the critical point 
the nucleon-nucleon and nucleon-non-Goldstone-meson 
interaction should gets much more attractive than it is in 
vacuum. We will now argue that it should affect collec- 
tive flows of nucleons, reducing their radial and elliptic 
flows. At the same time, we have argued in the previous 
chapter that the pion-pion interaction gets more repul- 
sive, so we expect the effect of the opposite sign in the 
pion flows. 

Potentials related to collective sigma field add to pres- 
sure effect in defining the flows. Vp term gets addition, 
rt s V<7. The distribution of the mean field sigma is to be 
found from the distribution of the scalar density via 
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(footnote: why not time derivatives? expansion is still 



slow) which should be solved together with hydro equa- 
tions. 

Treating its effect perturbatively 
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where the integrals are over the world line of each matter 
cell. 

Let us make a very simple order-of-magnitude estimate 
for the additional contribution to the pion flow we expect. 
The oval of fig.l projected to a fireball will correspond 
to also an ellipsoid-like surface. The pions crossing from 
inside out would experience a kick from the change in 
potential, repulsive inside and attractive outside. Non- 
relativistically, an extra velocity this kick will produce 
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where the final value in the r.h.s. is obtained by using 
the change in the potential A V e ff = 40MeV from Fig. 6 
and the effective pion mass for motion in one direction 
m 2 e ff =< p\ > +m\ ps (2T c ) 2 . This should be compared 
to collective radial velocity of matter at SPS, which also 
happen to be of the same magnitude, < v r >^ .5: the 
conclusion is the two effects are comparable. Thus ex- 
pected the non-monotonous change in pion radial flow 
can be quite noticeable. It will be also accompanied in 
increased elliptic flow of comparable magnitude. 

Similar estimate can be made for the nuclcon would 
include the potential change*** AU ~ .2GeV from Fig. 5 
and miff =< p% > +m 2 N w 1 GeV, resulting in 



Avpf ps —.5 — .6 



(17) 



Again, this is comparable to overall matter flow: but 
since the sign is now negative the total effect expected is 
cancellation of flows. Needless to say, any uncertainties 
in AU and schematic estimate like that cannot substitute 
for a detailed calculation: we how conclude that there is 
at least a potential for dramatic decrease of both radial 
and elliptic nucleon flows. 



A. Flows observed at SPS 

Let us now compare these ideas with observations. As 
it is well known, elliptic flow changes sign inside the AGS 
energy domain and is quite well seen at its highest energy, 
11 GeV/N. The nucleon and pion elliptic flows are also 
rather large at the highest SPS energy 158 GeV/N, and 
about double at RHIC. 



"'This is similar to what is done e.g. for magnetohydrody- 
namics, where one component of all fields is singled out from 
the rest of the matter, and its evolution is treated separately, 
coupled to the stress tensor and hydrodynamics. 



* M Note that the usual nuclear potentials would not be 
enough. 
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The summary of the 158 GeV/N elliptic flow data from 
NA49 experiment is shown in Fig. 7. The magnitude is 
enhanced at mid-rapidity, both for pions and nucleons. 
As is the case in other cases, and the pion and nucleon 
V2 are comparable in magnitude. 
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FIG. 8. Elliptic flow V2 of protons versus rapidity from 40 
A-GeV Pb+Pb collisions [16] measured for three centrality 
bins: central (dots) , mid-central (squares) and peripheral (tri- 
angles). The solid lines are polynomial fits to the data [16]. 

It is however not so at 40 GeV/N, as the same NA49 
experiment finds Fig. 8. Such collapse of flow at mid- 
rapidity is seen for all centralities and methods. It is 
furthermore accompanied by a collapse of directed flow 
Vi'. this coincidence tells us that it is unlikely to be an 
experimental problem and is probably real. 

Bratkovskaya et al. [27] suggested that this unusual 
behavior of the N flow may be a long-expected manifes- 
tation of the "softest point" . However that contradicts 
to the fact that elliptic flow of pions does not show any 
collapse at this energy. 

Our current proposal is instead that flows are affected 
by the approaching critical point. The nucleon flow col- 
lapse may be due to attraction effect we discussed, while 
pions show the opposite potential. 

Needless to say, this work is only exploratory in nature, 
and its suggestion should be investigated quantitatively 
in future work. In particular, to study flows, one should 
not only include space-time-dependent sigma field and 
potentials discussed above, but also cascades which in- 
clude continuous reactions like N + it— > A— > N + ir. 
Another issue worth considering is the effect of the criti- 
cal mode on strangeness flows, by <j>, A, S, VL. 

Finally, in view of all these ideas and observations, 
from theory, lattice and even hints from data, it seems 
justified to revisit the low SPS energy region. It can 
either be done with old SPS detectors (NA49,NA60) in 
new dedicated runs, or at decelerated RHIC beams, or 
maybe in future GSI dedicated facility to be built in the 
next decade. 
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